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Abstract. We write out the explicit form of the metric for a linearized 
gravitational wave in the transverse-traceless gauge for any multipole, thus 
generalizing the well-known quadrupole solution of Teukolsky. The solution 
is derived using the generalized Regge-Wheeler-Zerilli formalism developed by 
Sarbach and Tiglio. 



PACS numbers: 04.20.-q, 04.20.Jb, 04.25.-g, 04.25.D-, 04.30.-w 

1. Introduction 

In the hnearized approximation to general relativity, one assumes that the spacetime 
metric gap can be regarded as a small perturbation of the flat-space (Minkowski) 
metric rjap, 

where |(5ga/3| <^ 1. Solutions to the linearized Einstein equations have been known 
since the early days of general relativity. One usually works in the harmonic (or 
Lorentz, or de Bonder) gauge, in which the perturbation satisfies 

y^Sg^p - ^VaSg = 0. (1) 

Here V is the flat-space metric connection, Sg = Sg^p, and indices are raised and 
lowered with 77. Harmonic gauge does not fix the coordinates completely; one can 
impose four additional conditions. A popular choice is transverse-traceless (TT) gauge, 

Sgoc = 0, (2) 
6g = 0, (3) 

where denotes the time component. (Only four of these equations are independent 
of the harmonic gauge condition ([!]).) Conversely, if one only imposes ([2]) as a gauge 
condition then the components Roa = of the vacuum Einstein equations imply that 
the second time derivative of Q and the first time derivative of ([T]) also hold. 

A common way to represent the solutions is as expansions in tensor spherical 
harmonics, which are eigenfunctions of the (tensor) Laplace operator on the sphere 
(see [1] for a review article and [2] for an alternative approach). Teukolsky [3J wrote 
out the quadrupole {£ = 2) solution, giving all the functions in explicit form. This is 
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particularly helpful for numerical relativists who need expressions that can readily be 
coded. Teukolsky's has become one of the most widely used solutions, both for code 
testing purposes and for the construction of vacuum initial data using the conformal 
method (the two earliest references being [HIS]). As far as I know, such an explicit form 
of the solution has not appeared in the literature for ^ > 2. Having such a solution 
at hand will be useful in order to model higher-multipole gravitational waves. For 
instance, it can be used as a testbed for improved absorbing boundary conditions [6]. 
Being able to evolve the higher multipoles correctly is important, e.g. the contribution 
of octupole {£ = 3) radiation is found to be substantial in non-equal-mass binary black 
hole mergers [7j |8] . 

In section [2l we present the general solution of the linearized vacuum Einstein 
equations in TT gauge. The radial functions are written out explicitly. The angular 
functions are given in terms of spin-weighted spherical harmonics or alternatively in 
terms of derivatives of the standard spherical harmonics. The number of independent 
angular functions needed is reduced considerably as compared with T. Explicit 
expressions in terms of elementary functions are listed for £ = 2,3,4 in |Appendix A[ In 
section[3l we describe how the solution was derived using the improved gauge- invariant 
Regge- Wheeler- Zerilli O [10] formalism developed by Sarbach and Tiglio llj, which 
proves to be a powerful method for generating solutions of the linearized Einstein 
equations. 

2. The solution for arbitrary 1^2 

There are two independent polarization states, or equivalently even- and odd-parity 
waves. The even-parity metric is 



The functions Y , Yg, Y^, Ygo and Yg^ depend on 9 and (j) only; they are defined below 
in (fTO| or equivalently (fTTj) . The odd-parity metric is 

g = -df + dr^ + 2KSg r dr d6' + 2X84, r sin 61 dr d^ + (1 + LSggY"^ dO'^ 



Again the functions K and L depend on r and t only and are defined in ([9]). The 
functions Sg,S^, Sgg and S'e^ are directly related to the even-parity angular functions 
via (|12|) . Note that when written in this form, the perturbation is manifestly traceless, 
and we have reduced the number of independent angular functions to five as opposed 
to (seemingly) twelve in 

The solutions are specified in terms of mode functions F{x) for even parity and 
G{x) for odd parity that can be chosen arbitrarily. We set F^''^{x) = d'^ /dx''F{x) and 
similarly for G. Outgoing (ingoing) solutions are obtained by taking the argument of 
the functions F''"'-' and G''"'-' to be r — t (r + t). (By using r — t instead of t — r as in 
[3] we achieve the same form of the radial functions ([5]) and © for both directions of 
propagation.) We define the coefficients cj , ^ j ^ £, recursively by 




+2LSgy sin6ld6'd0+ (1 - LSggy"^ iiin^ dcp^ . 



(5) 



= 1 



J(2^-J + 1) 
2(£ - J + 1) 



Cj for j = 1, 



(6) 



Explicit solution of the linearized Einstein equations in TT gauge for all multipoles 3 
or equivalently, 

The radial functions appearing in the even-parity metric ([4]) are now given by 

e 

A = J2 CjT^^'^'^i-^i^ - jyF^^+^^ + 4(3/ - je^ - 7j£ + j + e^ + - 2£)rF(J+i) 
i=o 

+ (2/ - Ajl + e^- £){2j 2)f(J')], 

I 

B^J2 Cjr^^'^-^mi - j)r2F(^+2) - {6f - 12jl - 2j + 51^ + £ - 2)rF^'+^'> (8) 

3=0 

-{2f~Aje + e^-£)ij-e^2)F^% 

C = Y^ Cjr^'~^"^h2r2F(j'+2) - 2(2j - 2£ + l)ri^(j'+i' - {2f- - ijt + - e)F'^^'>], 

and those in the odd-parity metric (O are 

e 

3=0 

e 

L = J2 Cjr^"^~2[2j,G(J+i) + 2{j - £ + 1)G^^'>]. (9) 

3=0 

For a muhipole £ we need the derivatives of the mode functions up to order £ + 2 for 
even parity and up to order £ + 1 for odd parity. 

The angular functions in ^ are most easily written in terms of spin-weighted 
harmonics [12] syem{S,(j))- These can be obtained recursively from the standard 
spherical harmonics Yfm by 

oYem — Yim, 

s+iYem = [{£ ~ s){£ + s + sYe^, 

3_iy,„. ^~[{£ + s){£ - s + i)]-i/2e 

where the operators 9 and 9 are defined by 

dsYijn = {-de -icsc9d^ + s cot 9)sYera, 

dsYijn = {-de + icsc9d^- s cot 9)sYira- 
For simplicity we omit the indices £ and m. We have 

Y = oY, 

Ye = -ye(i+l){iY - -lY), 



Y4,^^^/£{£TT}(,Y + ^iY), (10) 

Yee = 1^(£-1)£{£+1){£ + 2){2Y + ^^Y), 

Ye^ = -i^{£-m£+m + 2){2Y - ^2Y). 

If preferred the following expressions in terms of the standard spherical harmonics and 
their partial derivatives may be used, 
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Y = Y, 

%^csc0Y^, (11) 
Yee = Yee + W + l)Y, 
Ye^ = CSC 9{Yfi^ - cot 6* F^). 
The odd-parity angular functions in ([5]) are related to the even-parity ones by 

Sb — —Ycf,, 5*0 — Yg, See = —Yg^, Sg^ — Ygg. (12) 

For each i there are 2i+l independent real solutions. These are obtained by replacing 
Yijn with KeYifn ii m ^ and with ImY^i^i if m < 0. The mode functions F and G 
must then be taken to be real. 



3. Derivation of the solution 

We have derived the solution using the generalized Regge-Wheeler-Zerilli l9l [10] 
formalism of Sarbach and Tiglio [11| . In this formalism, the background spacetime 
is assumed to be a direct product M ® S"^ of a Lorentzian 2-manifold M and the 
2-sphere. The metric on M is denoted by its volume element by e and its associated 
covariant derivative by V. In our case g is flat, 

g = -dt^ +dr2, 

and V reduces to the partial derivative. The standard metric on S"^ is denoted by g, 
its volume element by e and its associated covariant derivative by V. We have 

g = d6'2 + sin^ 61 d0^ 

The metric perturbation (in an arbitrary gauge for the time being) is decomposed 

as 

5gab = HabY, 

SgAb = QbYA + gbSA, (13) 
SgAB = r\KgABY + GYab) + 2kSAB. 

Here lower-case Latin indices a,b, . . . range over t and r and upper-case Latin indices 
A, B, . . . range over 9 and 0. The even-parity basis tensor spherical harmonics are 
derived from the standard scalar spherical harmonics Y according to 

Ya = VaY, 

Yab = [VaVbF]" = VaVbY + ^£{£ + l)gABY, 
where tf denotes the tracefree part. The odd-parity basis harmonics are 

Sa = aYb, 
Sab = ^{aSb) = [Sab]^^- 
Indices £ and m have been omitted and a sum over these in (jl3|) is implied. 
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3.1. Master equation 

Solutions of the linearized Einstein equations can be described in terms of two gauge 
invariant scalars 'f'^^^ the (generalized) Regge- Wheeler- Zerilli (RWZ) scalars. Here 
+ refers to even and — to odd parity. They obey the master equations 

[-V'^Va + l^±(r)]$W. 

For a flat-space background, V±{r) = £{£ + l)/r^ and hence 

iie + iy 



d - dt 



= 0, (14) 



an example of the Eulcr-Poisson-Darboux equation [13] . This equation is easily solved 
by making the ansatz 

i 

$(±)(t,r) ^^Cjr''^F^l\x) (15) 
i=o 

where F^{x) = /dx^ F±{x) and the argument is taken to be a; — r — t for an 
outgoing solution and x — r + t for an ingoing solution. The F± are precisely the 
mode functions appearing in the final solution (section [2]), with F = F+ and G = F-. 
For both directions of propagation (fTi|) implies 

ji2e - J + i)c, + 2{e - J + i)c,_i = 0, 

and hence together with the convention = 1 we obtain ^ and ([7]). 

We need to reconstruct the metric perturbation from the RWZ scalars and impose 
the TT gauge. The two parities are treated separately. Throughout we use the 
notation ' = d/dt and ' = d/dr. The reader is referred to ll for further explanation 
of the formalism underlying the following calculation. 



3.2. Odd parity 

From the Regge- Wheeler scalar <J>'^~\ we first compute the gauge-invariant potential 



(inv) 



eafcV"(r$ 



The bracket denotes the components [h^^'^^\ /ir™^-*] of the one- form hl'"\ The gauge- 
invariant potential is related to the amplitudes of the perturbation via 



r^Vb 



The gauge condition ^ implies ht — and so we obtain first 

k ^ 

and finally 



kV 



We define the integral to simply lower the index of the F^ by one, and in the case 
of the outgoing solution multiply the result by —1. No integration constant is being 
added. 
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3. 3. Even parity 

From the Zerilli scalar we obtain the ZeriUi one-form 



Zn 



AV„$ 



{+) 



where A = — + 2). Next we form the gauge-invariant potential 

r X r X 

From this the one-form C is defined by 



Ca=Za+ rVaK 



(inv) 



The gauge-invariant potential -ff^™^' can be deduced using 



and the fact that -ff^b"^'' tracefree. We find 



(inv) 



o-(inv) _ rj 



(inv) 



The metric amplitude Hab is related to this via 

where are gauge parameters that we compute first. The gauge condition ([2]) implies 
Htt = Htr = and hence subsequently 

2 



Pt 



Hir'dt, 



Pr = 
TT _ rrCinv) 



{Hfr'+p'^)dt, 



2^. 



The gauge parameters pa can also be expressed in terms of the amplitudes of the 
perturbation as 



Qt 



Sr-'G, Qr - 



The gauge condition ^ implies Qt = and hence subsequently 
2 



G 



Pr dt, 



Pr 



Finally 



K = + 2p^VV - U{t + l)G = xC'^^) + ^pr - U(i + 1)G. 
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3.4^. Final form of the solution 

Starting from the solution (jlSp of the master equation and going through the 
calculation described above, we arrive at the radial functions (O and ([5]), where we 
identify 

A^Hrr, B=^, C = G, K=^, L=^. 

Because of the integrations involved we have raised the derivative index of the F^^^ by 
2 and of the G^^'^ by 1 so that no negative derivative indices appear in the final form 
of the solution. The modified angular functions are defined by 

Y4, = CSC OY^, 84, = CSC e S^, Ye^ = esc 6 Yg^, 894, = esc 9 S04, 
and the remaining ones are equal to their unhatted counterparts. 

In our derivation of the solution, we have only imposed the gauge condition 
It can be checked directly that the remaining conditions ([T]) and ([3]) are also 
satisfied. This also follows more generally from the argument given below equation 
([3]) in the introduction, noting that our definition of the time integral merely lowers 
the derivative index of the while leaving the expressions formally unchanged. 
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Appendix A. Explicit expressions for ^ = 2, 3,4 



The (real) angular functions are listed below in the order m — £,£ — 1, . . . , 1, 0. We 
use the shorthand 



Rn 



cos mq 
sin md 



— sm mq 
cos md) 



where the two elements of each of these correspond to the two independent real 
solutions for each m > (the real and imaginary parts of the spherical harmonics). 
The odd-parity angular functions are immediately obtained from the even-parity ones 
using (112). For simplicity we leave out the normalization factors of the spherical 
harmonics. 



For I ~ 2 we have 

A =2A 

B = 4 
C = 2 
K = A 



F(2) 3^(1 



^3 

r 

G(2) 



L = 2 



^ 3£CT 

^ 2F(3) 
3G(i) 
2G(2) 



3F 

r* 
3f(^) 

3G 

3G(i) 



6F" 

3^ 

r4 



3F 



3G 
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Y ^sm^eR2, cos6lsin6li?i, 

Yg = 2 cose sin0R2, (1 - 2 sin^ 6')i?i, -6 cos 61 sin 6*, 

Y^ = 2sm9l2, cos9Ii, 0, 

Yge = {2 — sin'^ 6) R2, — cos 6* sin 6'i?i , 3sin^0, 

Ye4, = 2cos6'/2, -sin6l/i, 0. 

This solution agrees with that of Teukolsky [3] up to an overall constant factor 
(depending on m and the parity). 



For £ — 3 we have 

A = 120 

B =10 

C = 2 



F(3) Qp(2) 15^(1) i^p 

— ~ + ~ — ^ 

6i^(3) 21F(2) 45F(i) 45F 



5F(4) 15F(3) 30F(2) 45i^(i) 45F' 



X = 10 
L = 2 



G(3) 6G(2) 15G(i) 15G 



G(4) 5G(3) 15G(2) 30G(i) 30G 



r 



,.2 



y = sin 6'i?3, cos 6' sin 6'i?2, sin 0{4:-5sm'^0)Ri, cos6'(2 - 5sin^6'), 
Fe = 3cos6'sin2 6'i?3, sin6'(2 - Ssin^ 6')i?2, cos 61(4 - 15 sin^ e')i?i, 

3sin6'(5sin^6i-4), 
Y^ = 3sm'^ Oh, 2cos6lsin6l/2, (4 - 5 sin^ 6l)/i, 0, 

Fee = 3sin6'(2-sin2 6i)i?3, cos 0(2 - 3 sin^ 6l)i?2, 5sin6l(3sin2 - 2)i?i, 
15 cos 6* sin 6*, 

Fe^ = 6cos6'sin6'/3, 2(1 - 2 sin^ 6')/2, -10cos6lsin6'/i, 0. 



A = 360 
B =18 
C = 2 
K = 18, 
L = 2 



Finally, for ^ = 4 we have 

10F(3) 45i^(2) 105F(i) 105i^ 



- F^5) 


10i^(4) 


55F(3) 


195F(2) 


420F(i) 


420F' 

r' 












1 






45i^(4) 


150F(3) 


360i^(^ 


) 

- + 


630i^(i) 


630i^" 


r 
















"GW 


10G(3) 


45G(2) 


105G(i) 


105G' 












^4 




j,6 


1 







£-(5) 9G(4) 45G(3) 150G(2) 315G(i) 315G 

^ + ^ + ; ^ 
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F = sm''6'i?4, cos6isin^6ii?3, sin^ 61(6 - 7 sin^ 6')i?2, cos6'sin6'(4 - Tsin^ 6l)i?i, 

35 sin"* 61-40 sin^ 61 + 8, 
Ye = 4cos6isin^ OR^, sin^ 6'(3 - 4sin^ 6i)i?3, 4 cos 9 sin 61(3- 7 sin^ 6')i?2, 

(28 sin^ 6* - 29 sin^ 6i + 4)i?i , 20 cos sin 61(7 sin^ 61 - 4) , 
%=Asm^eh, Scosesin^eis, 2 sin6'(6 - 7 sin^ 6')/2, cos6'(4 - 7sin2 e*)/!, 0, 
Fee = 6sin2 6l(2-sin2 6l)i?4, 6 cos^ 6* sin 0i?3 , 6(7sin4 - 8sin2 6* + 2)i?2, 

6cos6lsin6'(7sin2 9 - 3)i?i, 30sin2 0(6 - 7sm^ 9), 
Ye^ = 12coaesm^9h, 3 sin 6'(2 - 3 sin^ 6l)/3, 6 cos 6'(2 - 7 sin^ 0)/2, 

3sin6'(7sin2 6'-6)/i, 0. 
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